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Abstract 

A model-independent search for deviations from the Standard Model prediction is per- 
formed using the full e^p data sample collected by the HI experiment at HERA. All event 
topologies involving isolated electrons, photons, muons, neutrinos and jets with transverse 
momenta above 20 GeV are investigated in a single analysis. Events are assigned to ex- 
clusive classes according to their final state. A dedicated algorithm is used to search for 
deviations from the Standard Model in the distributions of the scalar sum of transverse 
momenta or the invariant mass of final state particles and to quantify their significance. 
Variables related to angular distributions and energy sharing between final state particles 
are also introduced to study the final state topologies. No significant deviation from the 
Standard Model expectation is observed in the phase space covered by this analysis. 
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1 Introduction 



At HERA electrons^ and protons collide at a centre-of-mass energy of up to 319 GeV. The 
collected luminosity of high-energy electron-proton interactions gives access to rare processes 
with cross sections of the order of 0.1 pb, providing a testing ground for the Standard Model 
(SM) complementary to e + e~ and pp scattering. 

A large variety of possible extensions to the SM predicts new phenomena which may appear 
at high energies. Searches for new physics often compare the data to the predictions of specific 
models. A complementary approach is followed in signature based searches by looking for 
differences between data and SM expectation in various event topologies. As an advantage, 
such model independent analyses do not rely on any a priori definition of expected signatures 
for exotic phenomena. Therefore, they address the important question of whether unexpected 
phenomena may occur through a new pattern, not predicted by existing models. Following this 
approach, final states corresponding to rare SM processes such as single W boson or lepton 
pair production have already been investigated at HERA [1-5]. Model independent analyses 
are also performed at the Tevatron [6,7]. 

The present paper reports on a general analysis of all high transverse momentum (Pt) final 
state configurations involving electrons (e), muons (/i), jets (j), photons (7) or neutrinos (v) 
in e ± p collisions. This analysis searches for deviations from the SM prediction in phase space 
regions where the SM prediction is reliable. All final states containing at least two particles^ 
(e, /i, j, 7, u) with P T > 20 GeV in the polar angled range 10° < 9 < 140° are investigated. 
The present analysis follows the strategy of the previous HI publication [8]. Selected events are 
classified into exclusive event classes according to the number and types of particles detected 
in the final state (e.g. e-j, [i-v-j, j-j-j-j)- In a first step the event yields are compared with 
the SM expectation. In a second step kinematical distributions are systematically investigated 
using a dedicated algorithm [8] which locates the region with the largest deviation of the data 
from the SM prediction. 

The complete e ± p data sample collected by the HI experiment at HERA is used. The data 
are recorded at an electron beam energy of 27.6 GeV and proton beam energies of 820 GeV or 
920 GeV, corresponding to centre-of-mass energies ^/s of 301 GeV or 319 GeV, respectively. 
The total integrated luminosity of the data is 463 pb -1 , which represents a factor of four increase 
with respect to the previously published result [8]. The data comprise 178 pb -1 recorded in e~p 
collisions and 285 pb -1 in e + p collisions, of which 35 pb -1 were recorded at y/s = 301 GeV. 
While the previous general search was dominated by e + p collision data, a large data set recorded 
in e~p scattering is now also analysed. 

1 In this paper the term "electron" is used generically to refer to both electrons and positrons, unless otherwise 
stated. 

2 In this context a high Pt jet is also called particle. 

3 The origin of the HI coordinate system is the nominal ep interaction point, with the direction of the proton 
beam defining the positive z-axis (forward region). The transverse momenta are measured in the xy plane. The 
pseudorapidity -q is related to the polar angle 9 by r\ = — In tan(6*/2). 
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2 Standard Model Processes and their Simulation 



A precise estimate of all processes relevant at high transverse momentum in ep interactions is 
needed to ensure a reliable comparison to the SM. Several Monte Carlo (MC) generators are 
therefore combined to simulate events in all classes. 

At high transverse momenta the dominant SM processes are photoproduction of two jets 
and neutral current (NC) deep-inelastic scattering (DIS). Direct and resolved photoproduction 
of jets and prompt photon production are simulated using the PYTHIA [9] event generator. The 
simulation is based on Born level hard scattering matrix elements with radiative QED correc- 
tions. The RAPGAP [10] event generator, which implements the Born, QCD Compton and 
boson gluon fusion matrix elements, is used to model NC DIS events. QED radiative effects 
arising from real photon emission from both the incoming and outgoing electrons are simulated 
using the HERACLES [11] program. In RAPGAP and PYTHIA, jet production from higher 
order QCD radiation is simulated using leading logarithmic parton showers. Hadronisation is 
modelled with Lund string fragmentation [9] . The leading order MC prediction of photopro- 
duction and NC DIS processes with two or more high transverse momentum jets is scaled by 
a factor of 1.2 to account for the incomplete description of higher orders in the MC genera- 
tors [8, 12]. Charged current (CC) DIS events are simulated using the DJANGO [13] event 
generator, which includes first order leptonic QED radiative corrections based on HERACLES. 
The production of two or more jets in DJANGO is accounted for using the colour-dipole- 
model [14]. Contributions from elastic and quasi-elastic QED Compton scattering are simu- 
lated with the WABGEN [15] generator. Contributions arising from the production of single 
W bosons and multi-lepton events are modelled using the EPVEC [16] and GRAPE [17] event 
generators, respectively. 

All processes are generated with at least ten times the integrated luminosity of the data 
sample. Generated events are passed through the GEANT [18] based simulation of the HI 
apparatus, which takes into account the running conditions of the different data taking periods, 
and are reconstructed and analysed using the same program chain as is used for the data. 

3 Experimental Conditions 

A detailed description of the HI experiment can be found in [19]. Only the detector components 
relevant to the present analysis are briefly described here. The Liquid Argon (LAr) calorime- 
ter [20] covers the polar angle range 4° < 9 < 154° with full azimuthal acceptance. Electromag- 
netic shower energies are measured with a precision of a(E)/E ~ 11%/v/E/GeV © 1% and 
hadronic energies with a(E)/E ~ 50%/a/-E/GcV©2%, as measured in test beams [21,22]. In 
the backward region, energy measurements are provided by a lead/scintillating-fibre (SpaCal) 
calorimeter [23] covering the range 155° < 9 < 178°. The central (20° < 9 < 160°) and 
forward (7° < 9 < 25°) inner tracking detectors are used to measure charged particle trajec- 
tories and to reconstruct the interaction vertex. The innermost central proportional chamber, 
CIP [24,25] (9° < 9 < 171°) is used together with tracking detectors to veto charged particles 
for the identification of photons. The LAr calorimeter and inner tracking detectors are en- 
closed in a super-conducting magnetic coil with a field strength of 1.16 T. From the curvature 
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of charged particle trajectories in the magnetic field, the central tracking system provides trans- 
verse momentum measurements with a resolution of o> T / ' P T = 0.005Pr/GeV © 0.015 [26]. 
The return yoke of the magnetic coil is the outermost part of the detector and is equipped with 
streamer tubes forming the central muon detector (4° < 9 < 171°). In the forward region 
of the detector (3° < 9 < 17°) a set of drift chambers detects muons and measures their 
momenta using an iron toroidal magnet. The luminosity is determined from the rate of the 
Bethe-Heitler process ep — > ejry, measured using a photon detector located close to the beam 
pipe at z = —103 m, in the backward direction. 

The main trigger for events with high transverse momentum is provided by the LAr calorime- 
ter [27]. Events with an electromagnetic deposit (electron or photon) in the LAr with an energy 
greater than 10 GeV are detected by the LAr trigger with an efficiency of about 100% [28]. 
Events are also triggered by jets only, with a trigger efficiency above 95% for Pf et > 20 GeV 
and nearly 100% for Pf et > 25 GeV [29]. For events with missing transverse energy of 20 GeV, 
the trigger efficiency is about 90% and increases above 95% for missing transverse energy 
above 30 GeV [30]. The trigger for events with only muons is based on single muon signatures 
from the central muon detector, combined with signals from the central tracking detector. The 
trigger efficiency is about 95% for di-muon events with muon transverse momenta larger than 
15 GeV [5]. 



4 Data Analysis 

4.1 Event Reconstruction and Particle Identification 

In order to remove background events induced by cosmic showers and other non-ep sources, 
the event vertex is required to be within 35 cm in z of the nominal interaction point. In addition, 
topological filters and timing vetoes are applied [31]. 

Calorimetric energy deposits and tracks are used to look for electron, photon and muon 
candidates. Electron and photon candidates are characterised by compact and isolated elec- 
tromagnetic showers in the LAr calorimeter. The identification of muon candidates is based 
on a track measured in the inner tracking systems associated with signals in the muon detec- 
tors [1]. Calorimeter energy deposits and tracks not previously identified as electron, photon 
or muon candidates are used to form combined cluster-track objects, from which the hadronic 
final state is reconstructed [29, 32]. Jet candidates with a minimum transverse momentum of 
2.5 GeV are reconstructed from these combined cluster-track objects using an inclusive k T 
algorithm [33, 34] with a P T weighted recombination scheme in which the jets are treated as 
massless. The missing transverse momentum P™ 1SS of the event is derived from all detected 
particles and energy deposits in the event. In events with large P™ 1SS , a neutrino candidate is 
reconstructed. The four-vector of this neutrino candidate is calculated assuming transverse mo- 
mentum conservation and the relation ^(E* ~ p l) + ( E " ~ P z) = 2E e = 55 - 2 GeV > where 
the sum runs over all detected particles, P z is the momentum along the proton beam axis and 
P° is the electron beam energy. The latter relation holds if no significant losses are present in 
the electron beam direction. 
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Additional requirements are applied to ensure an unambiguous identification of particles, 
while retaining good efficiencies. Strict isolation criteria are applied in order to achieve high 
purities in all event classes. 

For electrons, the calorimetric energy measured within a distance in the pseudorapidity- 
azimuth (rj, 0) plane R = a/A?? 2 + A0 2 < 0.75 around the candidate is required to be below 
2.5% of its energy. In the region of angular overlap between the LAr and the central tracking 
detectors (20° < 9 < 140°), hereafter referred to as the central region, the calorimetric electron 
identification is complemented by tracking information. In this region it is required that a 
well measured track geometrically matches the centre-of-gravity of the electromagnetic cluster 
within a distance of closest approach (DCA) of 12 cm. Furthermore, the distance from the first 
measured track point in the central drift chambers to the beam axis is required to be below 
30 cm in order to reject photons that convert late in the central tracker material. In the central 
region, the transverse momentum of the associated electron track P^ tk is required to match 
the calorimetric measurement P| such that 1/ k — 1/Pf. < 0.02 GeV 1 in order to reject 
hadronic showers. In the forward region (10° < 9 < 20°), a wider calorimetric isolation cone 
of R < 1 is required to reduce the contribution of fake electrons from hadrons. In this forward 
region, at least one track is required to be present with a DCA < 12 cm. The presence of at least 
one hit in the CIP, associated to the electron trajectory, is also required. Finally, the electron 
is required to be isolated from any other well measured track by a distance R > 0.5 (R > 1) 
to the electron direction in the central (forward) region. The resulting electron identification 
efficiency is ~ 80% in the central region and ~ 40% in the forward region, determined from 
NC DIS events. 

The identification of photons relies on the same calorimetric isolation criteria as used in 
the electron identification. Vetoes on any track pointing to the electromagnetic cluster are ap- 
plied. No track with a DCA to the cluster below 24 cm or within R < 0.5 should be present. 
An additional veto on any hits in the CIP associated to the electromagnetic cluster is applied. 
Furthermore, each photon must be isolated from jets by R > 0.5. The resulting photon identi- 
fication efficiency as derived using elastic QED Compton events is ~ 95% in the central region 
and ~ 50% in the forward region. 

A muon should have no more than 5 GeV deposited in a cylinder, centred on the muon 
track direction, of radius 25 cm and 50 cm in the electromagnetic and hadronic sections of the 
LAr calorimeter, respectively. Misidentified hadrons are strongly suppressed by requiring that 
the muon be separated from the closest jet and from any further track by R > 1. In di-muon 
events, the opening angle between the two muons is required to be smaller than 165°, in order 
to remove muons originating from cosmic rays. The efficiency to identify muons is ~ 90% [5]. 

The scattered electron may be misidentified as a hadron and reconstructed as a jet. To reject 
fake jet candidates, the first radial moment of the jet transverse energy [35, 36] is required to be 
greater than 0.02 and the quantity M jct /P^ ct greater than 0.1 [12,36], where the invariant mass 
M jet is obtained using the four-vector sum of all particles belonging to the jet. If the fraction of 
the jet energy contained in the electromagnetic part of the LAr calorimeter is greater than 0.9, 
the above criteria are tightened to 0.04 and 0.15, respectively. These requirements are fulfilled 
by ~ 97% of the jets [36]. 

Missing transverse momentum, which is the main signature for neutrinos, may arise from 
mis-measurement of particles. By requiring £V (E l — P l z ) < 48 GeV, fake neutrino candidates 
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from NC DIS processes are rejected. If exactly one electron or muon candidate is found, a 
neutrino is only assigned to an event if A(f>n-x h ) < 160°, where A(f)n_x h ) is the difference in 
azimuthal angle between the lepton I and the direction of the hadronic final state X h . 

4.2 Event Selection and Classification 

The common phase space for electrons, photons, muons and jets is defined by 10° < 6 < 140° 
and P T > 20 GeV. The neutrino phase space is defined as missing transverse momentum above 
20 GeV and £\ (E { - P Z J < 48 GeV. All particles with P T > 20 GeV, including the neutrino 
defined by its reconstructed four-vector, are required to be isolated with respect to each other 
by a minimum distance R > 1. The particles satisfying these requirements are referred to 
as bodies. The events are sorted depending on the number and types of bodies into exclusive 
event classes. All possible event classes with at least two bodies are investigated. Only the 
fi-v event class is discarded from the analysis. This class is dominated by events in which a 
poorly reconstructed muon gives rise to missing transverse momentum, which fakes the neutrino 
signature. 

Based on these identification criteria, purities have been derived for each event class. Purity 
is defined as the ratio of SM events reconstructed in the event class in which they are generated 
to the total number of reconstructed events in this class. Most purities are found to be above 
60% and are close to 100% for the j-j, e-j, u-j and event classes. 

4.3 Systematic Uncertainties 

The following experimental systematic uncertainties are considered: 

• The uncertainty on the electromagnetic energy scale varies depending on the polar angle 
from 0.7% in the central region to 2% in the forward region. The polar angle measure- 
ment uncertainty of electromagnetic clusters is 3 mrad. The identification efficiency of 
electrons (photons) is known with an uncertainty of 3% (5%) to 5% (10%), depending on 
the polar angle. 

• The scale uncertainty on the transverse momentum of high P T muons is 2.5% [5]. The 
uncertainty on the reconstruction of the muon polar angle is 3 mrad. The identification 
efficiency of muons is known with an uncertainty of 5%. 

• The jet energy scale is known within 2% [30]. The uncertainty on the jet polar angle 
determination is 10 mrad. 

• The uncertainty on the trigger efficiency is estimated to be 6% if only muons are present 
in the final state and 3% in all other cases. 

• The luminosity measurement has an uncertainty of 3%. 
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The effects of the above uncertainties on the SM expectation are determined by varying 
the experimental quantities by ±1 standard deviation in the MC samples and propagating these 
variations through the whole analysis chain. 

Additional model uncertainties are attributed to the SM Monte Carlo generators described 
in section [21 An error of 10% is attributed to NC and CC DIS processes with only one high 
Pt jet. To account for the uncertainty on higher order QCD corrections, an error of 15% on 
the normalisation of NC DIS and photoproduction processes with at least two high P T jets is 
considered. The normalisation uncertainty of CC DIS processes with at least two high P T jets 
is estimated to be 20% [30]. For each additional jet produced by parton shower processes, a 
further theoretical error of 20% is added [37], for example 20% for the j-j-j event class. 

The error on the elastic and quasi-elastic QED Compton cross sections is conservatively 
estimated to be 5%. The error on the inelastic QED Compton cross section is 10%. The errors 
attributed to lepton-pair and W production are 3% and 15%, respectively. An uncertainty of 
30% on the simulation of radiative CC DIS events is considered to account for the lack of 
QED radiation from the quark line in the DJANGO generator. This uncertainty is estimated 
for the specific phase space of the analysis by a comparison of the DJANGO result to the 
calculated cross section of the e~p— ^u e ^X process [38]. An uncertainty of 50% is added to the 
prediction for NC DIS events with measured missing transverse momentum above 20 GeV and 
a high Pt electron. This uncertainty is estimated by a comparison of the missing transverse 
momentum distribution of data events containing a low P T electron (Pj, < 20 GeV) with the 
SM prediction [37]. 

The total error on the SM prediction is determined by adding the effect of all model and 
experimental systematic uncertainties in quadrature. 



5 Results 

5.1 Event Yields 

The event yields for all event classes are presented for the data and SM expectation in fig- 
ures \V[a) and UXb) for e + p and e~p collisions, respectively. All event classes with observed 
data events or with a SM expectation greater than 0.01 events are shown. The corresponding 
observed and predicted event yields for all e ± p data are summarised in table [TJ Events are ob- 
served in 27 classes and a good description of the number of observed data events by the SM 
prediction is seen in each class. 

The j-j, j-j-j and j- j-j-j event classes are dominated by photoproduction processes. No 
event with five jets is observed. The SM prediction of the e-j, e-j-j, e- j-j-j and e- j-j- j-j event 
classes is dominated by NC DIS processes. One event, already discussed in a previous HI 
publication [8], is observed in the e-j-j-j-j event class and compares to a SM prediction of 
0.13 ± 0.06. The u-j, v-j-j, v- j-j-j and v-j-j-j-j event classes mainly contain events from CC 
DIS processes. One event is observed in the v-j-j-j-j event class compared to a SM expectation 
of 0.05 ±0.02. 
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Events from QED Compton processes populate the 7-e event class as well as the 7-e-j event 
class in the case of inelastic events. The 7-j event class corresponds to prompt photon events. 
The purity in this class is moderate (~ 50%) due to the high background from misidentified 
electrons in NC DIS. A slight deficit of data events is observed in the radiative CC DIS classes 
7-1/ and 'y-u-j. 

Lepton pair production from 77 processes dominates in event classes with several leptons 
(e-e, fjL-fji, e-fj, and e-e-e). Compared to the results of a previous study of multi-lepton topolo- 
gies [5], the phase space of the present analysis is restricted to higher P T and extended to 
forward polar angles down to 10°. All multi-lepton events mentioned in [5] and located in the 
phase space of this analysis are found. The e-e event class contains 7 events with an invariant 
mass M ee > 100 GeV compared to a SM expectation of 3.4 ± 0.5 of which 69% are from 
lepton pair processes. The e-e-e event class contains one event compared to a SM expectation 
of 0.22 ±0.04. 

The prediction for the event classes [L-v-j and e-v-j consists mainly of high P T single W 
production with subsequent leptonic decay. In the [i-v-j (e-v-j) event class 5 (4) events are 
observed, with a SM expectation of 2.8 ± 0.5 (3.2 ± 0.5). Two events classified as n-v-j in the 
previous analysis [8] now migrate to fx-j and u-j event classes, respectively, due to improve- 
ments in the energy and momentum reconstruction. Events arising from W production also 
enter in the e-v event class. In this class 16 events are observed compared to an expectation of 
21.5 ± 3.5, of which about 90% is due to W production processes. 



5.2 Event Topology 

The distributions of the scalar sum of transverse momenta ^ Pt an d of the invariant mass M all 
of all bodies are presented in figures [2] and [31 respectively, for classes with at least one event. 
The data are in agreement with the SM prediction. In particular, multiple jets topologies, which 
are sensitive to QCD radiation, are well described by the simulation. 

The final state topologies are also evaluated in terms of angular distributions and energy 
ratios, which are sensitive to spin and decay properties of hypothetical high mass particles. 
Variables used to study the decomposition of the final states, inspired by topological analyses 
of multi-jet events [39], are defined in the following. In each event a leading body is selected 
according to the following priority list between bodies of different types: 7, e, //, v, j. This 
order of preference allows a better separation of SM background from events originating from 
a new resonance decaying to a photon or a lepton. If two bodies of the same type are present, 
the one with the highest transverse momentum P T , relative to the incident proton in the centre- 
of-mass frame defined by all bodies, is selected. For classes with exactly two bodies of the 
same type, the leading body is taken as the one with the highest Pt in the laboratory frame. The 
variable cos #j* ead is then defined as the cosine of the polar angle of the leading body relative to 
the incident proton in the centre-of-mass frame defined by all bodies. The variable Xi cad is the 
energy fraction of the leading body and is defined for systems with three or more bodies as 

9 W* 

V _ load r 1 \ 

-Mead - UJ 
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where the sum runs over all bodies energies, and E* ead and E* are calculated in the centre- 
of-mass frame of all bodies. For events with two bodies, the cos0* ead distribution is related to 
the underlying 2^2 matrix element. Therefore, the angular distribution of a particle coming 
from the decay of a new resonance may be markedly different from that of particles produced 
in SM processes (see for example [40]). For final states with more than two bodies, Xi ca d is a 
Dalitz variable and related to the dynamics of a possible multi-body decay of a new particle. 
The sensitivity of these two variables cos 9^ cad and Xi ead to new physics is tested using different 
MC samples of exotic processes, for example leptoquarks, excited fermions, or anomalous top 
production. It has been verified that SM and exotic events exhibit different spectra in these two 
variables, two examples of which are given in figure @] 

The distributions of cos0* ead and Xi ea d are presented in figure |5] for event classes with only 
two bodies and for event classes with more than two bodies, respectively. A good overall agree- 
ment with the SM prediction is observed in all cases. This illustrates that the event topology 
and kinematics, as well as the global variables Y Pt and M a n, are well described by the SM. 

5.3 Search for Deviations from the Standard Model 

In order to quantify the level of agreement between the data and the SM expectation and to 
identify regions of deviations in the P T , M all , cos0f ead and Xi ead distributions, the search 
algorithm developed in [8] is used. A region is defined as a set of connected histogram bins 
with at least twice the size of the resolution. A statistical estimator p is defined in order to judge 
which region is of largest interest. This estimator is derived from the convolution of the Poisson 
probability density function (pdf) to account for statistical errors and a Gaussian pdf to include 
the effect of systematic uncertainties [8]. The value of p gives an estimate of the probability 
of a fluctuation of the SM expectation upwards (downwards) to at least (at most) the observed 
number of data events in the region considered. The region of greatest deviation is the region 
having the smallest p-value, p m m- The regions selected by the algorithm in Y Pt and M a u 
distributions of each class are presented for all e ± p data in figures [2] and [31 respectively. The 
corresponding selected regions for cos 9* eBd and X load distributions are shown in figure [5J 

The fact that the deviation could have occurred at any point in the distribution is taken 
into account by calculating the probability P to observe a deviation with a p-value p min at 
any position in the distribution. P is a measure of the statistical significance of the deviation 
observed in the data. The event class of most interest in the search for anomalies is the one with 
the smallest P value. Values of P larger than 0.01 indicate event classes where no significant 
discrepancy between data and the SM expectation is observed. The P values measured in each 
of the event classes are listed in table \T\ Due to the uncertainties of the SM prediction in 
the j-j-j-j, e-j-j-j-j and u-j-j-j-j event classes, no reliable P values can be calculated for 
them [8] and they are therefore not considered in the search for deviations from the SM. 

The overall agreement with the SM can further be quantified by taking into account the 
large number of event classes in this analysis. Among all studied classes there is some chance 
that small P values occur. This probability can be calculated on a statistical basis with MC 
experiments. A MC experiment is defined as a set of hypothetical data histograms following 
the SM expectation with an integrated luminosity equal to the amount of data recorded. The 
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complete search algorithm and statistical analysis are applied to MC experiments analogously 
as to the data. The expectation for the P values observed in the data is then given by the 
distribution of P values obtained from all MC experiments. 

The P values observed in the data in all event classes are compared in figure |6] to the dis- 
tribution of P obtained from a large set of MC experiments. The comparison is presented for 
the scans of the M a n and Y Pt distributions for all e ± p data and also separately for e~p and 
e + p data. The distribution of P values measured in the data is in agreement with the expec- 
tation from MC experiments. Using all e ± p data, a lowest P value of 0.0044 is found in the 
e-j-j event class in a region at high transverse momenta, 175 < J^Pt < 200 GeV, where 27 
events are observed for an expectation of 11.6 ± 1.2. In e~p data, the lowest P value is 0.0071 
and corresponds to the e-e-e event class where one data event is observed compared to a low SM 
expectation. The most significant deviation from SM predictions is measured in e + p collisions 
in the e-e event class with P = 0.0035. In the corresponding region (110 < M a u < 120 GeV) 
five data events are found while 0.43 ± 0.04 are expected. The global probability to find in the 
e + p data at least one class with a P value smaller than observed in the e-e event class is 12% as 
deduced from MC experiments. 

In case of the cos 9* egd and Xi eaxi distributions, no significant discrepancy between the data 
and the SM expectation is found. The lowest P value is 0.017, observed in the Xi ead distribution 
of the v-j-j-j event class. In event classes where the SM contribution is high (> 100 events), 
the correlation between M all or Y P T distributions and cos 9* ead and X lead is further exploited. 
The variables cos 8* ead and X lea d are used to select events in a phase space region where the SM 
contribution is reduced and exotic event topologies may be favoured. Events where the leading 
body is emitted in the forward direction are selected by requiring cos 0f ead > 0. The variable 
Xie a d is used in three bodies event classes to select topologies corresponding to a sequential 
resonance decay by requiring 0.75 < Xi ca d < 0.9, as deduced from the study of different 
MC samples of exotic processes. After a cut on these variables, an overall good agreement 
between the data and the SM is still observed in M all and Y Pt distributions. The complete 
search procedure and statistical analysis is applied to these distributions, the results of which 
are summarised in tabled No significant deviation is observed in the reduced event samples. 

The full analysis is also performed at lower and higher transverse momenta by changing the 
minimum P T of particles to Pt > 15 GeV and Pt > 40 GeV, respectively. A good overall 
agreement with the SM is also observed with these cuts. With a cut Pt > 15 GeV, all spectra 
are well described by the MC, including the multi-jet event classes. The lowest P value is 0.01, 
observed in the e-j-j event class. When raising the Pt threshold to 40 GeV, mainly event classes 
containing jets remain populated and the largest deviation is observed in the e-j-j-j class with 
P = 0.01. 



6 Conclusion 

The full e ± p data sample collected by the HI experiment at HERA is investigated in a general 
search for deviations from the SM prediction at high transverse momenta. This analysis en- 
compasses all event topologies involving isolated electrons, photons, muons, neutrinos and jets 
with transverse momenta above 20 GeV. Data events are found in 27 different final states and 
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events with up to five high P T particles are observed. In each event class deviations from the 
SM are searched for in the invariant mass and sum of transverse momenta distributions using 
a dedicated algorithm. In addition, the final state topologies are also evaluated in terms of an- 
gular distributions and energy sharing between final state particles. A good agreement with the 
SM expectation is observed in the phase space covered by this analysis. The largest deviation 
is found in the e-e event class, in e + p collisions, at high invariant masses and corresponds to 
a probability of 0.0035. The probability to observe a SM fluctuation with that significance or 
higher for at least one event class is 12%. This comprehensive analysis demonstrates the very 
good understanding of high P T SM phenomena achieved at the HERA collider. 
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Table 1: Observed and predicted event yields for all event classes with observed data events 
or a SM expectation greater than 0.01 for all e ± p data. Each event class is labeled with the 
leading body listed first. The errors on the predictions include model uncertainties and experi- 
mental systematic errors added in quadrature. The P values obtained in the scan of P T , M all , 
cos 9* ead and X lcad distributions are also given. 
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HI General Search at HERA (e±p, 463 pb" 1 ) 
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Table 2: Observed and predicted event yields for considered event classes after a cut on the 
topological variables. Each event class is labeled with the leading body listed first. The errors 
on the predictions include model uncertainties and experimental systematic errors added in 
quadrature. The P values obtained in the scan of Pt and M a n distributions are indicated in 
the last two columns. 
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H1 General Search at HERA (e + p, 285 pb" 1 ) H1 General Search at HERA (e"p, 178 pb" 1 ) 
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Figure 1: The data and the SM expectation for all event classes with observed data events or a 
SM expectation greater than 0.01 events for e + p collisions (a) and e~p collisions (b). The error 
bands on the predictions include model uncertainties and experimental systematic errors added 
in quadrature. 
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H1 General Search at HERA (e^p, 463 pb 1 ) -EP T Distributions 
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Figure 2: The number of data events and the SM expectation as a function of Pt f° r classes 
with at least one event. The shaded areas show the regions of largest deviation identified by the 
search algorithm. No such search is performed for the j-j-j-j, e-j-j-j-j and u-j-j-j-j classes. 
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H1 General Search at HERA (e^, 463 pb" 1 ) - M aM Distributions 





100 200 300 100 200 300 100 200 300 



M all [GeV] 

Figure 3: The number of data events and the SM expectation as a function of M all for classes 
with at least one event. The shaded areas show the regions of largest deviation identified by the 
search algorithm. No such search is performed for the j-j-j-j, e-j-j-j-j and u-j-j-j-j classes. 
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Figure 4: The cos6* cad distribution in the 7-e event class (a) and the Xi cad distribution in the 
e-j-j event class (b). The points correspond to the observed data events and the histograms 
to the SM expectation. The error bands on the SM prediction include model uncertainties 
and experimental systematic errors added in quadrature. The dashed line represents, with an 
arbitrary normalisation, the distribution corresponding to an exotic resonance with a mass of 
200 GeV (e* [41] in (a) and v* [42] in (b)). 
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H1 General Search at HERA {erp, 463 pb" 1 ) 
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Figure 5: The distribution of cos for event classes with two bodies (top) and of X lead for 
event classes with more than two bodies (bottom). The points correspond to the observed data 
events and the open histograms to the SM expectation. Only event classes with at least one data 
event are presented. The shaded areas show the regions of largest deviation identified by the 
search algorithm. No such search is performed for the j-j-j-j, e-j-j-j-j and v-j-j-j-j classes. 
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H1 General Search at HERA (e ± p, 463 pb" 1 ) 
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H1 General Search at HERA (e'p, 178 pb" 1 ) 
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H1 General Search at HERA (e ± p, 463 pb" 1 ) 
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Figure 6: The -log 10 P values for the data event classes and the expected distribution from MC 
experiments as derived with the search algorithm by investigating the M a u distributions (left 
column) and Pt distributions (right column). The results of the scan is presented for all data 
(a and b), and separately for e'p (c and d) and e + p (e and f) data. 
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